UNCLASSIFIED _ 

itciwrv  CLASSIFICATION  OF  This  PAGE 


DOCUMENTATION  PAGE 

AD-A234  082 

lb  RESTRICTIVE  MARKINGS 

3  .  DISTRIBUTION  /AVAILABILITY  OF  REPORT 

Approved  for  public  release;  distribution 
is  unlimited. 

2b.  DECLASSIFICATION /DOWNGRADING  SCHEDULE 

4.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

USAFSAMhJA*90*27 

5.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

6*.  NAME  OF  PERFORMING  ORGANIZATION 

6b  OFFICE  SYMBOL 
(If  applicable) 

7a.  NAME  OF  MONilORING  ORGANIZATION 

nTi.fi . 

6c  ADDRESS  (City,  State,  and  ZIP  Code) 

7b.  ADDRESS  (City,  State,  and  ZIP.CP#)  ,  .  .  .  ^  .  .. 

8a.  NAME  OF  FUNDING  /  SPONSORING 
ORGANIZATION 

Armstrong  Laboratory 

8b.  OFFICE  SYMBOL 
(If  applicable) 

Det  4  AL/RZM 

9  PROCUREMENT  INSTRUMENT  IDENTIFICATION, DUMBER  t|gj| 

Sc  ADDRESS  (City.  State,  and  ZIP  Code) 

Armstrong  Laboratory  (AFSC) 

Brooks  AFB  TX  78235*5301 

10.  SOURCE  OF  FUNDING  NUMBERS 

PROGRAM 
ELEMENT  NO. 

62202F 

PROJECT 

NO. 

2729 

TASK 

NO 

5 

WORK  UNIT 
ACCESSION  NC 

03 

1 1 .  TITLE  (Include  Security  Classification) 

AnimaluTo?1Huraan  Extrapolation  Using  Compartmental  Models 


12.  PERSONAL  AUTHOR($) 

Medina,  Richard  L.;  Albanese,  Richard  A. 


13a.  TYPE  OF  REPORT  j 

13b.  TIME  COVERED 

14.  DATE  OF  REPORT  (Year,  Month,  Day) 

15  PAGE  COUNT 

Interim  1 

FROM  10/89  TO 

04/90 

1990,  April 

18 

16.  supplementary  notation 


J  17.  COSATI  CODES  | 

£  FIELD 

GROUP 

SUB-GROUP 

rTT 

I  1 

_ 

04  1 

jf 

I  r.-Air* 

J  v'S. 

i  -***» 
I  * 

i 


-I 


18.  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  u y  block  number) 

Fickian  Equation,  Venous  Exit  Condition,  Compartmental 

Modeling,  Diffusion  Model,  Mass  Balance,  Extrapolation 
Partial  Differential  Equation 


19.  ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

We  have  been  studying  how  to  use  compartmental  models  to  reliably  extrapolate  toxicologic; 
experiment  results  from  one  animal  species  to  another  with  the  ultimate  goal  of  useful 
extrapolation  to  man.  We  have  taken  as  the  fundamental  core  of  our  analysis  the  physical: 
necessary  equations  describing  mass  balance.  These  equations  are  the  classical  Fickian 
equations  of  the  form  VdC/dt  =  Q(a-v).  These  equations  are  not  a  complete  set  since  they 
are  insufficient  to  provide  estimation  of  venous  arterial  or  tissue  concentrations  with 
time.  Therefore,  the’ mass  balance  equations  must  be  augmented  with  additional  phenomenon 
logical  relations  to  permit  the  desired  calculations.  Many  investigators  use  the-  venous 
exit  condition.  While  it  is  clear  how  to  extrapolate  the  mass  balance  equations  from  one 
sized  animal  to  another  or  from  one  species  to  another,  it  is  not  clear  how  to  extrapolat 
the  venous  exit  condition.  In  this  research  we  have  evaluated  the  venous  exit  condition 
by  comparing  it  with  approximate  analyses  of  perfusion  and  substance  diffusion  in  the 
tissue. 


JO.  DISTRIBUTION/  AVAILABILITY  OF  ABSTRACT 

K3  UNCLASSIFIED/UNLIMITED  □  SAME  AS  RPT.  Q  DTlC  USERS 


21 


ABSTRACT  SECURITY  CLASSIFICATION 

Unclassified 


J2».  NAME  OF  RESPONSIBLE  INDIVIDUAL 

Medina,  Richard  L. _ 


22b  TELEPHONE  ( Include  Are*  Code) 

51 2P‘536«-3S8h 


22c  OFFICE  SYMBOL 
Det  4  AL/RZM 


DD  Form  1473,  JUN  86 


Previous  editions  are  obsolete 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 

- tmCLfiSSIPTED - 


Neuroscience  <4  Biobehaviorat  Reviews.  Vol.  15.  ^p,  57-61.  Pcrgamon  Press  pic.  1991.  Primal  in  the  U.S.A. 


01 49-76  U/91  S3 .00  +  .00 


Animal-to-Human  Extrapolation 
Using  Compartmental  Models 


RICHARD  L.  MEDINA  AND  RICHARD  A.  ALBANESE 

United  Stales  Air  Force  School  of  Aerospace  Medicine 
Brooks  AFB,  San  Antonio,  TX  78235-5301 


MEDINA,  R.  L,  AND  R.  A.  ALBANESE.  Animal-to-liuman  extrapolation  using  compartmental  models.  NE'IROSCI  B10BEHAV 
REV  15(1)  57-61.  1991.  — We  have  been  studying  how  to  use  compartmental  models  to  reliably  extrapolate  toxicological  experi¬ 
mental  results  from  one  animal  species  to  another  with  the  ultimate  goal  of  useful  extrapolation  to  man.  We  have  taken,  as  the  fun¬ 
damental  core  of  our  analysis,  the  physically  necessary  equations  describing  mass  balance.  These  equations  are  the  classical  Fickian 
equations  of  the  form  VdC/dt  -  Q(a  -  v).  These  equations  are  not  a  complete  set  since  they  are  insufficient  to  provide  estimation 
of  venous,  arterial,  or  tissue  concentrations  with  time.  Therefore,  the  mass  balance  equations  must  be  augmented  with  additional 
phenomenological  relations  to  permit  the  desired  calculations.  Many  investigators  use  the  venous  exit  condition.  While  it  is  clear 
how  to  extrapolate  the  mass  balance  equations  from  one  sized  animal  to  another  or  from  one  species  to  another,  it  is  not  clear  how 
to  extrapolate  the  venous  exit  condition.  In  this  research  we  have  evaluated  the  venous  exit  condition  by  comparing  it  with  approx¬ 
imate  analyses  of  perfusion  and  substance  diffusion  in  the  tissue. 

Fickian  equation  Venous  exit  condition  Compartmental  modeling  Diffusion  model  Mass  balance 
Extrapolation  Partial  differential  equation 


THIS  report  describes  research  the  authors  have  pursued  concern¬ 
ing  compartmental  modeling  in  toxicology  with  particular  regard 
to  the  problem  of  animul-to-mari  scale  up.  Several  benefits  would 
accrue  if  medical  scientists  could  reliably  extrapolate  biological 
effects  findings  from  one  sized  animal  to  another,  or  from  one 
species  to  another.  For  example,  imagine  that,  in  mice,  a  chemi¬ 
cal  administered  for  3  months  at  a  daily  dost  D  results  in  an  in¬ 
crement  in  cancer  C,  where  C  **  f(D),  when  the  presence  or 
absence  of  cancer  is  determined  by  autopsy  performed  when  the 
mice  reached  2  years  of  age.  Given  the  data  based  equation,  C  = 
l‘(D).  for  mice,  what  relation  between  cancer  and  exposure  might 
exist  for  man?  The  question  concerning  rnan  would  usually  go 
beyond  a  regimented  daily  dose  D  given  only  for  a  standard  3- 
month  period,  and  concern  for  cancer  would  extend  for  the  entire 
human  life  span. 

This  article  addresses  the  introduced  class  of  problems.  The 
approach  applied  to  this  problem  is  compartmental  modeling.  We 
believe  the  prospects  for  this  endeavor  can  be  favorable  as  is  al¬ 
ready  indicated  in  the  pharmacokinetics  literature.  In  this  article, 
we  will  address  the  general  methodology  of  extrapolation  using 
compartmental  modeling,  and  will  not  be  treating  a  specific  ap¬ 
plication. 

TWO  PLAUSIBLE  APPROACHES  TO  COMPARTMENTAL  MODELING 

For  purposes  of  discussion  we  will  use  the  block  diagram 
shown  in  Fig.  1 ,  Figure  1  represents  a  simplified  version  of  Ram¬ 
sey  and  Andersen's  PB-PK  model,  which  described  the  pharma¬ 
cokinetics  of  styrene  in  rats  and  humans  during  and  after  inhalation 
exposure  (9). 

In  the  setting  of  this  block  model,  the  following  mass  balance 
equations  are  applicable: 

VL  dCL/dt  =  QL(v  -  a) 


v„  dC„/dt  =  Q„(a  -  v„) 

VM,  dCuJ/dt  =  Qad(a  -  vad) 
dCjhs/dt  ~  Qabd3  —  vabsl  hahs(\h  —  Cabs) 

V.xc  dCc,,/dt  =  Qe,c(a  -  vCXL.)  -  \C„, 

Qv  =  Qnv„  +  Qad  vad  +  Qabs  vab,  +  Q„c  v„t 

(EQUATION  1) 

The  five  symbols,  Vj,  represent  compartment  volumes  and  the 
five  symbols,  Qj,  are  compartment  flows.  The  symbols  C,  refer 
to  substance  concentration  in  each  organ.  The  following  subscript 
convention  is  used:  L  designates  the  lung,  n  designates  nonadi- 
pose  tissue,  ad  designates  adipose  tissue,  abs  designates  the  ab¬ 
sorbing  organ  (for  example,  the  intestine  or  integument),  and  exc 
designates  the  excretory  organ  (for  example,  the  kidney),  The 
symbol  a  refers  to  the  arterial  concentration  of  substance,  v  refers 
to  the  combined  venous  concentration  and  the  symbols  v,  refer  to 
concentrations  of  material  in  the  venous  outflow  of  the  individ¬ 
ual  organs  (excluding  the  lung).  The  symbol  Xs  corresponds  to  an 
external  substance  concentration  impinging  on  the  absorbing  sur¬ 
face. 

These  equations  are  all  a  direct  reflection  of  the  principle  of 
conservation  of  matter  (3).  Except  for  the  kabs(xs  -  Cab,t  and 
\Ccxl.  terms,  which  are  not  necessarily  applicable,  these  equa¬ 
tions  must  be  true  in  any  compartmental  system  because  they  are 
based  on  physical  principles.  In  these  equations  the  V,,  Q,,  kab5 
and  \  terms  are  considered  known,  but  the  Cr  a.  v  and  v,  terms 
tire  not.  Thus,  there  are  eleven  unknown  variables  and  only  six 
equations.  Therefore,  the  question  .irises  of  how  to  complete  these 
equations,  or  close  the  set.  so  the  Cr  v  and  Vj  values  can  be  cal- 
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FIG.  1.  A  five  compartment  block  diagram  of  a  physiological  pharmaco¬ 
kinetic  model  studied  in  this  paper. 


culated  given  Vj*  Qj'  katn  and  X. 

One  way  to  complete  the  mass  balance  equation  set  is  through 
use  of  the  venous  exit  condition  (8).  This  condition  asserts  that 
blood  flowing  out  of  an  organ  always  has  a  substance  concentra¬ 
tion  proportional  to  the  concentration  of  the  substance  in  the  or¬ 
gan.  The  following  venous  exit  condition  equations  can  be  written 
to  complete  the  mass  balance  set. 


a  '  aLCL 


an  Cn 

«e*c  Cexv 


A 


»  f  t 


X=0 


FIG.  2.  The  proposed  perfusion -diffusion  mode!  showing  substance  dif¬ 
fusion  into  a  tissue  segment  supplied  by  a  central  blood  vessel. 


^7ad  dCjKj/cit  —  kad  [C^  (a  +  vad)/2] 
v„bl  dCab,/dt  =  -kabs  icabs  -  (a  +  vabJ/2]  +  k,b,  (x5  -  C,ba) 
V„c  dC„c/dt  =  -k„t.  [C„c  -  (a  +  v„v.)/2]  -  *C„« 


(EQUATION  3) 

Thus,  we  now  have  two  complete  compartmental  models  cor¬ 
responding  to  Fig.  1 .  A  part  of  each  model  is  physically  exact, 
being  based  on  conservation  of  matter.  However,  each  model  also 
contains  mathematical  expressions  that  do  not  exactly  reflect 
physical  laws,  but  are  approximations.  Specifically,  the  approxi¬ 
mations  in  the  models  use  algebraic  or  ordinary  differential  equa¬ 
tions  when  the  perfusion-diffusion  phenomena  in  tissue  is  actually 
governed  by  partial  differential  equations  (8).  Two  questions  im¬ 
mediately  arise.  First,  which  of  the  two-compartmental  models 
will  best  describe  substance  pharmacokinetics  in  a  living  animal? 
Second,  which  of  these  two  models  is  most  accurate  when  used 
to  extrapolate  from  a  smaller  to  a  larger  animal? 

The  best  way  to  answer  these  two  questions  is  through  the  use 
of  extensive  data  obtained  from  different  sized  animals,  for  ex¬ 
ample.  rats  and  humans.  Since  the  large  data  sets  needed  were 
not  available  to  us  we  have  pursued  an  alternative  analysis  of  the 
competing  compartmental  models.  We  have  compared  the  venous 
exit  and  diffusion  approximation  models  to  a  simple  perfusion- 
diffusion  system  governed  by  partial  differential  equations  for 
which  analytic  mathematical  solutions  exist.  This  test  perfusion- 
diffusion  system  exhibits,  qualitatively  at  least,  much  of  the  ac¬ 
tual  substance  movement  behavior  of  tissue.  We  have  determined 
how  well  the  compartmental  models  approximate  this  exact  per¬ 
fusion-diffusion  system,  and  we  have  studied  how  well  the  com¬ 
partmental  models  follow  a  scale-up  of  the  exact  system. 


(EQUATION  2) 

The  coefficients  a,  may  be  thought  of  as  simply  constants  of 
proportionality  or  may  be  thought  of,  more  specifically,  as  in¬ 
versely  proportional  to  the  partition  coefficients  (9).  The  venous 
exit  condition  is  a  simple  algebraic  expression  which  may  or  may 
not  hold  true  in  any  given  case  This  condition  is  an  attempt  to 
approximate  a  complex  perfusion-diffusion  system  actually  gov¬ 
erned  bv  partial  differential  equations  (8). 

It  is  possible  to  complete  or  close  the  mass  balance  equations 
in  a  second  way.  Recognizing  that  transport  into  an  organ  com¬ 
partment  occurs  through  a  diffusion  process,  approximate  diffu¬ 
sion  equations  may  be  written  (6,7).  We  have  used  the  following 
diffusion  approximation  closure  equations: 


AN  EXACT  PERFUSION-DIFFUSION  SYSTEM 

We  propose  to  examine  and  develop  an  exact  solution  of  the 
Sangren  and  Sheppard  model  of  tissue  clearance  (4,10),  We  con¬ 
sider  a  Krogh  cylinder  (4,  10,  11)  of  length  L  as  diagrammed  in 
Fig.  2.  The  tissue  segment  is  supplied  by  a  central  blood  vessel 
with  volume  -y  and  the  concentration  of  the  substance  in  the  blood 
B(x,t),  varies  as  a  function  of  linear  position  and  time.  Let  C(x,t) 
be  the  concentration  in  the  tissue  itself.  Assuming  simple  diffu¬ 
sion  across  the  vessel  wall  and  no  axial  or  radial  diffusion  in  the 
tissue,  the  following  equations  apply: 

•ydB(x.t)/tit  =  k[C(x.t)  -  B(x.t)]  -  -yp,dB(x,t)/dx 
VdCtx.tj/Tt  =  -kjC(x.t)  -  B(x.l)] 


VL  dC,  /dt  ~  -kL  (CL  -  (a  +  v)/2]  (EQUATION  4) 

Vn  dCn/dt  =  -kn  [Cn  -  (a  +  vn)/z]  In  these  equations,  B  =  B(x.t)  is  the  concentration  of  substance 


AN1MAL-TO-HUMAN  EXTRAPOLATION 


5*) 


B(0,t) 


H 


{l£|(4:>  +  id  +  i,  (4i  +  fell  sin(t*Jjt) 

+  [?i  (“t*  ~  id  +  i,  (£2  +  iO]  COS(ujjt)} 


►  t 

H/2  D  2D-H/2  2D 

FIG.  3.  The  substance  temporal  profile  in  the  blood  entering  the  Kmgh 
cylinder. 


in  the  blood,  C  =  C(x,t)  is  the  concentration  of  substance  in  the 
tissue.  7  is  the  volume  of  blood  in  the  tissue,  p  is  the  velocity  of 
blood  flow  in  the  tissue,  V  is  the  tissue  volume,  and  k  is  the  per¬ 
meability  constant  for  the  vessel  wall. 

We  have  taken  B(o,t)  to  be  the  trapezoidal  function  of  time  t 
shown  in  Fig.  3.  This  function  has  the  following  Fourier  series 
representation  1 1 ): 


X 


B(o,t)  =  (1/2)  -  4/[tt"(1- 


(1/j3)  cos(j-rrp)  cos  (u^t) 


where  =  jir/D  and  p  -  H/2D.  Using  this  Fourier  series  rep¬ 
resentation.  equation  4  can  be  readily  solved  to  find 


B(x,t)  =  (1/2)  -  4/[-tr3(l  -2p)]  Y 
(l/j3)  exp(-t|ijX/p)  cos(jirp)  co$(u>j(t - 9,)] 


and 


C(x.t)  =  (1/2)  -  4(k/V)/[Tr-(  1  -  2p)] 


where 


€.  =  o>/[«V  +  MV)2] 

£.2  =  l  I/(Otj~  +  (Uljpj)-)]  [-Oj  COSflUjpL) 

+  (lOjpjtsinluJiBjLllexpf-ajL) 

£3  ~  «/[« j2  +  Ibl.Pjl'l 

i»  =  1 1/(0, 3  +  (WjBj)3))  1-Oj  sinfto.PjL)  -  (uijPj) 
cosfuijpjLijexpl  -  OjL) 

ks  =  (lOjPjt/lQj3  +  (to^j)3] 

6,  -  (k/v)/[tOj3  +  (L'Y)3; 

“j  =  ttylA 

pj  =  (1/p)  {to/  +  (k3/V)  [(1/V)  +  (l/7)]}/(ojj3  +  (k/V)3) 

In  the  exact  perfusion-diffusion  model,  ave[C(t)]  is  the  value  cor¬ 
responding  to  the  Cj  terms  in  the  compartmental  approximations. 
Also  B!o.t)  corresponds  to  substance  concentration  in  in-flowing 
blood,  while  B(L,t)  corresponds  to  the  concentration  in  tissue 
exiting  blood. 

COMPARING  THE  EXACT  PERFUSION  DIFFUSION  SYSTEM  TO  THE 
VENOUS  EXIT  AND  APPROXIMATE  DIFFUSION  COMPARTMENTAL 
MODELS 

For  a  specific  compartment,  the  venous  exit  model  is 
VdC/dt  +  Q  a  C  =  0  a 
v  =  a  C 


(l/j3)  exp  (-ill,  X/p)  cos(jnp). 

{|u>J/(u>J2  +  (k/V)2)]  sin|u>,(t  -  Bjl]  +  |th/V)/(tu3  +  (k/Vl2)] 
cos(u>j(t  —  8,)]} 

where 

^  =  w,3  (k/yl/Ou,3  +  (k/V)3) 

Gj  =  (x/p)  {u>j3  +  (k3/V)[( i/V )  +  (l/ylU/tajj3  +  (k/V)3) 

Mean  tissue  substance  concentration  ave|C(t)]  can  be  calculated 
as: 


In  this  model,  wc  'ake  a  =  B(o.t)  and  solve  for  C  and  v  to  be 
compared  with  avelC(t)]  and  B(L.i)  respectively  from  the  exact 
perfusion-diffusion  system. 

Similarly,  for  a  specific  single  compartment,  the  approximate 
diffusion  model  is 

VdC/dt  =  Q(a  —  v) 

VdC/dt  =  -k[C  -  (a  +  v)/2] 

These  equations  are  equivalent  to 

VdC/dt  +  ]2kQ/(2Q  +  k)]C  =  [2kQ/(2Q  +  k)]a 


ave]C(tl)  =  (1/L)J^  Ctx.tldx 


v  =  [2k/(2Q+  k)]C  +  [( 2Q  -  k)/(2Q  +  k)]a 


whence 


avelC(t)]  =  (1/2)  -  4(k'rV)/(-iT2L(  1  - 


Again,  we  take  a  =  B(o.t)  and  solve  for  C  and  v  to  compare  with 
uve|C(t)]  and  B(L,t)  respectively  from  the  exact  perfusion-diffu- 

si  on  system. 

We  have  contrasted  the  performance  of  the  venous  exit  com¬ 
partmental  model  and  the  approximate  diffusion  compartment 
model  to  the  exact  perfusion-diffusion  system  by  treating  rat  and 
human  liver.  For  the  rat  liver  V  =  12  cm3  and  Q  =  35.2  cm’/ 


( !/j2)  cos(jtrp). 
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FIG.  4.  Fitted  data  showing  slow  change  in  rat  venous  blood  in  all  three  FIG.  6.  Extrapolation  results  for  man  regarding  the  estimation  of  venous 
models.  concentration. 


min.  For  the  human  liver  V  =  3320  cm’  and  Q  =  1800  cm’/min 
(9),  We  have  used  V  =  3.14  x  10  5  cm’  and  -y  =  7.85  x 
10“*  cm’  in  the  exact  perfusion-dif  usion  model,  assuming,  for 
test  purposes  only,  that  both  rat  and  human  liver  are  made  up  of 
cylindrical  regions  I -mm  long,  and  100  p.m  in  radius  with  a  cen¬ 
tral  blood  vessel  having  a  5  p.m  radius  (2,  4,  5,  1 1).  For  the  ex¬ 
act  model,  blood  velocity  was  chosen  to  be  117  cm/min  in  the  rat 
so  the  exact  model  would  have  a  Q/V  ratio  of  35.2/12,  and  blood 
velocity  was  chosen  to  be  21.7  cm/min  in  the  human  to  corre¬ 
spond  to  a  Q/V  ratio  of  1800/3320. 

Figures  4  through  7  show  the  results  of  our  numerical  experi¬ 
mentation.  In  these  figures  best  fit  a  and  k  values  were  found  so 
that  the  venous  exit  and  approximate  diffusion  compartments  best 
mimicked  the  rat  exact  model.  Extrapolating  the  human  liver  sit¬ 
uation,  a  was  left  unchanged  while  k  was  extrapolated  using  the 
relationship 

K2km  Q„,)/(2Qr„,  +  krnl)  VrJ  =  [(2kman  Q „,„„)/ 

(2Qma„  +  km„n)  Vmon] 


In  Fig.  4  we  see  that  all  three  models  can  be  made  to  agree  as 
regards  a  slow  change  in  rat  venous  blood.  However,  in  Fig.  5 
we  see  that  the  same  parameters  which  gave  agreement  in  rat 
venous  blood  do  not  work  as  well  for  rat  tissue  concentrations.  In 
particular,  the  venous  exit  condition  provides  a  poorer  quality  fit 
to  rat  tissue  concentrations. 

To  correspond  to  man.  the  exact  perfusion-diffusion  system 
blood  velocity  was  changed  from  117  cm/min  to  21.7  cm/min. 
The  parameter  a  was  left  unchanged  for  the  venous  exit  model 
and  k  was  adjusted  as  described  before.  In  Fig.  6  we  see  that  the 
venous  exit  model  and  diffusion  approximation  model  both  ex¬ 
trapolate  well  to  man  as  regards  the  estimation  of  venous  concen¬ 
trations.  However,  as  seen  in  Fig.  7,  the  venous  exit  condition 
again  does  not  perform  well  regarding  tissue  concentrations. 

The  patterns  found  in  Figs.  4  through  7  are  amplified  when 
fasier  changes  in  blood  levels  are  studied.  In  these  settings  again 
tile  diffusion  approximation  model  excels  the  venous  exit  model. 
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FIG.  5.  Rat  tissue  concentration  for  all  three  models.  Same  parameters  FIG.  7.  Extrapolation  results  for  man  regarding  the  estimation  of  tissue 
used  which  gave  agreement  in  rat  venous  blood.  concentration. 
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CONCLUSION 

Compartmental  models  have  a  core  set  of  mathematical  rela¬ 
tions  which  are  derived  from  the  physical  principle  of  conserva¬ 
tion  of  matter.  This  core  set  of  equations  is  not  sufficient  for 
modeling,  but  must  be  augmented  by  auxiliary  conditions.  The 
venous  exit  condition  is  commonly  used,  and  a  diffusion  approx¬ 
imation  also  is  possible.  These  approximations  use  algebraic  equa¬ 
tions  or  ordinary  differential  equations  to  describe  perfusion- 
diffusion  phenomena  :n  tissue  governed  by  partial  differential 
equations.  We  compared  a  venous  exit  model  and  a  diffusion  ap¬ 
proximation  model  to  an  exact  perfusion-diffusion  model  of  tis¬ 
sue  governed  by  partial  differential  equations.  We  have  found 


that  the  diffusion  approximation  model  performs  better  than  the 
venous  exit  model  in  mimicking  the  perfusion-diffusion  system, 
particularly  as  regards  the  prediction  of  tissue  concentrations.  We 
were  surprised  to  observe  that  a  model  that  fits  blood  concentra¬ 
tion  changes  well  may  not  fit  tissue  concentration  events. 

Only  more  experience  with  a  variety  of  substances  and  real 
world  animai-to-human  extrapolation  efforts  will  provide  secure 
guidance  concerning  choice  among  competing  compartmental  mod¬ 
els.  This  report  suggests  that  comparison  of  competing  compart¬ 
mental  models  to  perfusion-diffusion  models  governed  by  partial 
differential  equations  may  be  an  additional  important  guide  to 
compartmental  model  selection. 
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